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■ Abstract — Rotation-induced splitting of the otherwise degenerate photonic bands is pre- 

dicted for a two-dimensional photonic crystal made of evanescently coupled microcavities. The 
£NJ ■ symmetry-broken energy splitting is similar to the Zeeman splitting of atomic levels or elec- 

^ ■ tron's (hole's) magnetic moment sublevels in an external magnetic field. The orbital motion 

of photons in periodic photonic lattice of microcavities is shown to enhance significantly such 
Coriolis- Zeeman splitting as compared to a solitary microcavity [D.L. Boiko, Optics Express 
2, 397 (1998)]. The equation of motion suggests that nonstationary rotation induces quantum 
transitions between photonic states and, furthermore, that such transitions might generate high- 
frequency gravitational waves. 

1. INTRODUCTION 

The Sagnac effect in a rotating ring cavity, known also as the Coriolis-Zeeman effect for photons, 
emerges as a frequency splitting of counterpropagating waves [T]. Thus, for a ring cavity with optical 
path length of M wavelengths, the modal shift is Mf2 with Q being a rotation rate. Nowadays, 
i i ' the effect is used in commercial He-Ne ring laser gyros of large (M~10 6 ) cavity size. Significant 

efforts were made towards designing miniature-sized solid-state devices based on a high optical gain 
medium [2]. Theoretical investigations have been carried out about the impact of rotation on the 
whispering gallery modes of a microdisc microstructure [3]. Here, the frequency splitting scales 
with the closed optical path length, while the field polarization is either not relevant or assumed 
, to be parallel to the rotation axis fi. 

On the other hand, the Coriolis-Zeeman effect in a cylindrical microwave resonator rotating 
(f-^ , along the symmetry axis [4J or in an optical Fabry-Perot cavity rotating in the mirror plane is, 

at first sight, independent of the cavity size [5]. The frequency shift (S+M)£l of polarization and 
0^ , transverse modes is set by the spin (±1) and the azimuth mode index M. However, the higher the 

mode index M, the larger the size of the cavity needed to support such a mode. By virtue of the 
complexity of the mode discrimination at high M and because of the small frequency splitting of 
the polarization modes, this effect has not yet been verified in experimental measurements. 

Here, the Coriolis-Zeeman effect is considered in coupled microcavities arranged in a periodic 
two-dimensional (2D) lattice in the plane of rotation [Fig Ufa)]. On example of a square-symmetry 
lattice, the possibility of enhanced Coriolis-Zeeman splitting, corresponding to M~1000, is pre- 
dicted for the low- order photonic modes [£]. It is caused by the photon's orbital motion extended 
over the large number of lattice cells. The equation of motion, which is similar to the Hamiltonian 
for electrons and holes in magnetic field, suggests that nonstationary rotation induces quantum 
transitions between photonic states and, vice verse, that such transitions will generate a nonsta- 
tionary gravitational field. 

2. PHOTONIC CRYSTAL OF COUPLED MICROCAVITIES 

Arrays of evanescently coupled microcavities belong to a particular sub-class of 2D photonic crystal 
(PhC) structures encompassing photonic crystal fibers and arrays of microcavities. Matrices of 
vertical cavity surface emitting lasers (VCSELs) are an example of such 2D photonic crystals [7]. 
In such structures, only a small transversal component kj_ of the propagation vector k undergoes 
Bragg reflections in the plane of periodic lattice [FigUfa)]. The structures employ lattices of periods 
significantly exceeding the optical wavelength. They are typically realized by the mirror reflectivity 
patterning in a broad-area microcavity. As a model system for such PhC structures, a Fabry-Perot 
cavity with patterned mirror reflectivity is considered here. The cavity length is one wavelength, 
l z =X/n with n being the refractive index in the cavity. The reflectivity Ri(x,y) of the one cavity 
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Figure 1: (a): Schematic of the array of coupled microcavities. (b): Band structure of square-lattice PhC 
calculated from [black curves] and using 8x8 kp approximation (fT0|) [gray curves]. The inset shows 
the BZ and location of the high-symmetry triangle AZT . The parameters are A=960nm, n=3.53, A=4/im, 
and A</?=0.02. 

mirror (e.g., of the top mirror) is modulated in two directions parallel to its plane. The reflectivity 
pattern Ri = exp iip(x, y) consists of pixels with the relative phase shift ip(x,y)=Aip separated by 
a grid of (p(x,y)=0. The mirror is thus perfectly reflecting (\Ri \ =1) and the phase modulation 
pattern (f(x, y) defines the structure of the cavity modes. The period of the reflectivity pattern A 
is of a few micron pitch (A~^>l z ). The pattern is characterized by a fill factor FF that is the ratio 
between the area of the pixel and that of the unit cell. Like in typical VCSEL arrays, the phase 
contrast of the reflectivity pattern is small (|A</?| ~10 -2 ). The second mirror of uniform reflectivity 
(i?2=l) has no impact on the cavity modes. 

3. MODEL HAMILTONIAN 

The analysis is carried out here using an equivalent, unfolded cavity representation [8]. Multiple 
reflections at the cavity mirrors effectively translate the cavity into a structure that is periodic along 
the cavity axis (the z-axis). The unfolded PhC is thus three-dimensional and it can be analyzed 
in terms of a modal expansion on orthogonal plane waves (OPWs). The Coriolis-Zeeman effect in 
such photonic crystal is considered here using a frame of reference, which rotates together with the 
crystal. Such noninertial rest frame is characterized by a metric tensor g^, the off-diagonal space- 
time components go a of which are dependent on the angular rotation speed fi [9]. However, in the 
unfolded-cavity representation, g^ differs from the diagonal Minkowski tensor, even in the absence 
of rotation. Thus goa-T^O a ^ the subsequent mirror reflections since the reflection operator a=IC2 
includes rotation by n (about the z axis) followed by coordinate inversion [FigJTJa)]. Due to the 
patterned mirror reflectivity, the equivalent unfolded PhC is of periodically varying " noninertiality" 
in the xy plane [8]. In the approximation of the first-order terms Qr/c and A<p, the metric tensor 
gik has the following nonzero components: 

5oo = -5n=-522=-533=l J (1) 
9oa = — e al3j Q p x^-5 a3 -ip(x 1 ,x' 2 )y^5(x 3 -2jl z ), 
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where go a =9aO (a = 1,2,3), the space-time coordinates are indexed according to the intervals 
dx°=cdt, dx 1 =dx, dx 2 =dy and dx 3 =dz; twice repeated Greek indexes indicate summation. The 
first term in go a accounts for rotation of the coordinate system pQ. The second term accounts for 
the multiple cavity roundtrips along the z-axis and reflections at the cavity mirrors. The z-period 
of the unfolded crystal is thus 21 z . The metric tensor ([TJ is validated by inspecting the system 
Hamiltonian [Eq.([7])] for the case of (p=0 (rotating FP cavity [5]) or $7=0 (PhC in an inertial frame 

0). 

In the approximation ([I]), the coordinate space is Euclidean, with the metric tensor 70-/3= — 
9a/3 + 9oa9o/3 / 9oo being the Kronecker delta 5 a p. Proceeding in a standard manner [9], the covariant 
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Maxwell's equations with metric ([T]) are converted to the usual form in terms of noncovariant field 
vectors B, H, D and E that assume the constitutive equations [TO] 
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D=eE + H x g, B=/iH + gxE, 
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where z is the unit vector along z-axis direction and the components of the vector g are g a = — 

dQa/dOO- 

Maxwell' equations in photonic crystal ([2]) are solved here by separating fast oscillations in the 
z-axis direction and slow lateral field oscillations in the xy plane: [TO] 
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where n=y / e/I and Z=y^i/e are the refractive index and impedance in the cavity. The gauge 
transformation is introduced here through the operator 

£ =5 (l 1 92 \ 1 92 

ap a/3 y ^2 f)n r _ f) nr _ ] 9^,2 f)T-f)n?-o 

(4) 
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where the terms ~£r[/&^ are taken into account. Such separation of variables is valid in conditions 



of the paraxial approximation ( 
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<Cl) and of the low contrast of reflectivity 
in the x-y plane is a slowly-varying 



pattern (|A</?| <Cl). The two-component vector i/'( r ±) = ^^ 
function of coordinates. It defines the spatial patterns of the six electromagnetic field components 
(j3J) and it is considered here as the photonic state wave function. Its squared modulus |-0(rj_) | 2 
yields the intensity pattern of the main polarization component in ©. For 17=0, Eqs.([3|)-(j4]) are 
in agreement with the results obtained for the Gaussian beam 

In photonic crystals, the wave function ip(r±) is a Bloch wave propagating in the xy plane [8], 



(5) 



where ^ j cM u* q , k± u qk± d 2 r ± =5 q , q [H [TD]. The longitudinal part in © [the term e ik * z l+ ^ z) ] 
is also a Bloch function. Within the z-period of the lattice, it has a small modulation depth 
( I I ) 2Z = 2F I-i \Vqk\ dz^Aip, which is set by an effective phase shift a q ^ at each reflection of the 
patterned mirror. The general form of such periodic function r/ g k is 



l+?7 g k= exp 



9(z-2jl z )-l 



2L 



(6) 



where 0(z)= f* 8(C)dC is the unit step function. Note that rjg^z) is the odd function and 
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iO by virtue of the small contrast of the reflectivity pattern. 



By operating with e^ a p£p and £ a ^ [from @] on Maxwell' equations for the curl of E and H, 
substituting the gauge ([3]) and averaging over the ^-period of the lattice, the Maxwell' equations 
are converted into the same form of a Hamiltonian eigenproblem with respect to the photonic state 
wave function t/> 9 k( r ±) 
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<p{r±) r_L x p_L + hS z 

n z 



lf>qk = TtUqkfpqk, 



(7) 



where rriQ=nhk z /c and S z =izx is the spin operator that reads (S z ) a p=ie a 3p= (9 l ) in the basis 
of the two-component vector functions ip= J. Solutions of Eq.([7|) define the slowly- varying 
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components of photonic modes ([3]). The difference between the exact equation for (l+Tfak^gk 
and its z-period average [Eq.flTJ)] yields the equation for the periodic part of the fast longitudinal 
component: 



= «« 9 k (1 + Vqk) 



(8) 



where a q k= (V\jk| V 9 l^ok)- The solution of ([8]) is given by ([6]) provided that the eigen functions ipq^ 
of the Hamiltonian are known. 

In (|7J), the first term is due to the paraxial propagation along the z axis. The second and third 
terms are the in-plane kinetic energy and the periodic crystal potential, respectively. The last term 
in j7|) is a perturbation induced by the Coriolis force. For (p=0 and 17=0, within the accuracy of 
the time variable, Eq.([7]) is just the paraxial wave equation. For ip=0, Eq.© yields the effective 

k 2 

refractive index ^(fca+^F") °^ a circularly-polarized paraxial wave 

fir 

n cS = n + gr± , (9) 

Lun 

where t=^ defines the propagation direction, +/— sign is for the left/right handed polarization. 
Eq.([9|) agrees with previously reported expressions for the axial nonreciprocity [1] (second term) 
and circular birefringence [5] (third term) induced by the Coriolis force for photons. Finally, the 
unperturbed ($7=0) Hamiltonian has been verified by experimental measurements in VCSEL array 
PhC heterostructures [12j . These justify the approximation ([7]) of the Hamiltonian. 

4. RESULTS AND DISCUSSIONS 

Analytical similarities between the effective single-electron Hamiltonian in a semiconductor sub- 
jected to an external magnetic field and Eq.([7]) allow the correspondence between the periodic 
crystal potential V and phase pattern <p (V— >— ^x^), and the vector potentials A=^Hxr and 
g_i_=^f2xr (| A— >^rgj_). Photons in photonic crystal subjected to a nonpermanent gravitational 
field thus exhibit a behaviour similar to electrons (holes) in a magnetic field. Accordingly, the 
impact of rotation on the envelope function and periodic part of the photonic Bloch wave © is 
different. As in the case of electrons, |13] the components of velocity operator vj_=^[r_|_, H] do 

not commute (\v T ,v v ] = ^%^ z , where v i = — — ^ x 2 rj - ). However, the second-order terms 

in the Hamiltonian (|7|) [and, respectively, in ([1])] are needed to define whether the Landau-like 
quantization is possible for photonic envelope wave functions. In the rest of the Letter, the impact 
of rotation on the periodic part of Bloch functions is examined in detail using the example of a 
square-lattice PhC. 

FigH{b) shows the typical band structure of a square-lattice PhC, which is calculated along the 
high symmetry lines A-Z-T in the Brillouin zone (BZ), using the OPW expansion in unperturbed 
(17=0) Hamiltonian. By virtue of the square lattice symmetry, all states are degenerate by the pho- 
ton's spin (e.g., the doubly degenerate states T5 or T' 5 ). Angular rotation removes the degeneracy 
of such states and splits their energies on 2^ [Figj2](a), top panel]. However, there are states, 



Table 1: Basis functions (scalars) and photonic harmonics (vectors) of irreducible representations of the 
group of k at the T point of the BZ (C^ v point group) 

_Ti Tj x T 5 

Ti S T 5 SZ,Sy 

T 2 XY(X 2 -Y 2 ) T' 5 ' XY(X 2 -Y 2 )Z,XY(X 2 ~Y 2 )y 

T 3 X 2 -Y 2 T' 5 " (X 2 -Y 2 )Z,(X 2 ~Y 2 )y 

T 4 XY % XY±,XYy 

T 5 iX,iY Ti+T 2 + T 3 + T 4 : 

Tx ^(xst + Yy) 

T 2 -^(YZ-Xy) 
T 3 jfiXi-Yy) 
T 4 j=(YZ + Xy) 
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Figure 2: (a): Coriolis-Zeeman splitting of the T5 (T'5) (top panel) and T1-T4 (bottom panel) bands as a 
function of rotation rate f2. A = Afim and A<p = 10~ 4 . (b): Relative splitting At ^" s (left axis) and matrix 
element (Ti|r 2 |Ti)^ (right axis) as a function of the lattice contrast Aip. The lattice constant A is 4 (black 
curves) and 6)im (gray curves), n = 3.53. Other parameters are given in the caption of FigfT] 



like the degenerate states Ti, T2, T3 and T4, of the four-fold degeneracy, which is caused by the 
orbital symmetry of the Bloch functions. For such states, there is an important angular momentum 
contribution \Ml z L z to the energy shift. 

The Coriolis-Zeeman splitting of these states is analyzed here using the first-order kp expansion 
in the T point [k = (t, ^, k z )\ of the BZ. The expansion basis is deduced from the empty lattice 
test. Four scalar plane waves exp(±^±^), which form the first photonic band of empty lattice, 
originate from the nearest equivalent T points of reciprocal lattice and provide representation that is 
reducible under the C/± v point group (the symmetry group of k). Their symmetrized combinations of 
Ti, T4 and T5 representations are indicated in the second column of Table Q] with the capitals letters 

corresponding to the main term of Taylor expansion on parameter (e.g., \S) =^ cos ^ cos t^, 
\iX) =— sin ^ cos ^). The photon's spin transforms as the two-dimensional representation T5. 
Therefore, a reduction of the direct product Ti®T§ results in the eight symmetry adapted photonic 
harmonics of T1-T5 and T' 5 representations that constitute a suitable kp-expansion basis in the 
low-order photonic bands. These states are well separated energetically from the other states in 
the T point [Fig. 1(b)]. 

For a general state IV'ok) =e l]l±r± ^ q |Tj) at kj_ measured from the T point of the BZ, the 8x8 
kp Hamiltonian for the coefficients Cj is of the block-diagonal form 
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in the basis of functions ^(T' 5jc ±iT^) , ^(T^iT,) , ±|(T 3 ±iT 4 ) , ±j(T 5x ±iT 5y ) [JJ]. The upper (lower) 



sign refers to the top (bottom) 4x4 block. In (fTUl) . k±=k x ± ik y , P- 
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V2 
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is the 
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interband matrix element of p that defines the band mixing. The eigen solutions of (jlOp are plotted 
in Fig. 1(b). Good agreement with the band structure calculated from Eq.© justifies the first-order 
kp approximation (|10|) of the system Hamiltonian. 

In Eq. (110j) . the matrix Hq in (llOj) is a perturbation induced by rotation. It accounts for the 
Coriolis-Zeeman energy shift, which is of the opposite sign for the left (upper 4x4 block) and right 
(lower block) handed polarization states of the same orbital symmetry. The rotation lifts degeneracy 
between such spin states, producing the frequency splitting Acos = 2^, as indicated in Figl^a) 
[top panel] on example of the T5 (T' 5 ) states. In (fTUj) . the parameter M= M + +M_ accounts for 
the orbital part of wave functions. The orbital contribution Acj£=2^P to the Coriolis-Zeeman 
energy shift ±^hAajL±^hAujs in the T1-T4 states is shown in the bottom panel of Figj2][a) (Aojg 
is not visible at the scale of Auj£). The energy shift ±^HAujl [the term — Arxp in Eq. ([?])] is 
evaluated here using the relationship r mn =j^ r p ™f i ■ Note that the procedures to evaluate the 
matrix elements of KL Z 



"m E n — E„ 

=r _i_ x P_i_ in f ree space [15] and in periodic lattices [16] are different. Here, 



the /-sum rule 



Jcp 1 2 
m ml 
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implies that M± = =f|(- 



-1) and 



Auj l 



n 2 



m m 



Auj s 



n 



(11) 



For an array of square pixels defining the microcavities, M±- 



2nl,P 2 



sin 7ry FF 
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It follows that reducing the effective mass (via the lattice pitch A, fill factor FF and contrast A(p), 
one can enhance the Coriolis-Zeeman splitting and achieve M of more than 10 3 . [FigE^b) shows 
left axis]. The enhancement is caused by the weak localization of photonic 

^^- fi2(M - + ^inthe 



2P 2 



the ratio ^ 

wave functions to the lattice sites. The intraband matrix element ((/k|r 2 |(/k) = 
T1-T4 bands indicates that the photonic wave functions spread over a large PhC crystal domain 
[Figj2]Jb) right axis], and the frequency splitting is proportional to the characteristic size of this 
domain 

2vr, 



Auj l _ 

n 



2(rj 



n 2 A 



[1+sinc 7ta 



(12) 



Thus, as in the case of a ring cavity, the frequency splitting of optical modes increases with the 
characteristic modal size in the plane normal to the rotation axis. Finally note that in the bands 
of T5 symmetry, the intraband matrix elements of r 2 are nonzero as well. However, the orbital 
contribution to the Coriolis-Zeeman energy splitting vanishes, in accordance with the group theory 
selection rules. 



5. CONCLUSION 

With present experimental techniques, the predicted frequency splitting [Eq.(|12j). Figl^a)] can be 
validated by direct measurements. 

The analogy between Eq.([7|) and electron's (hole's) magnetic moment Hamiltonian suggests 
that a nonstationary rotation or gravitational field g(t) will induce quantum transitions between 
photonic states. This effect might be used for detection of high frequency gravitational waves. 
Furthermore, Eq.([7]) allows formally an inverse process to take place as well. On the basis of such 
formalism, one can expect that a superposition of nonstationary photonic states, which represent a 
system undergoing quantum transition, might serve as a source g(t) for high frequency gravitational 
waves. [T7] 
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